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Abstract 

Pulse broadening for optical solitons due to birefringence is investigated. We present an analytical solution which 
describes the propagation of solitons in birefringent optical fibers. The special solutions consist of a combination 
of purely solitonic terms propagating along the principal birefringence axes and soliton-soliton interaction terms. 
The solitonic part of the solutions indicates that the decay of initially localized pulses could be due to different 
propagation velocities along the birefringence axes. We show that the disintegration of solitonic pulses in birefringent 
optical fibers can be caused by two effects. The first effect is similar as in linear birefringence and is related to 
the unequal propagation velocities of the modes along the birefringence axes. The second effect is related to the 
nonlinear soliton-soliton interaction between the modes, which makes the solitonic pulse-shape blurred. 
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I. Introduction 



The demand for long-distance high-bandwidth data transfer in the field of optical telecommunications has lead to 
research about the use of optical envelope solitons as information carriers. Solitons are ideal for optical telecommu- 
nications because the dispersion of the optical fiber is exactly counterbalanced by the nonlinearity. As a result of 
this, the information carrier can maintain its pulse-shape over long distances. In an ideal optical fiber, the nonlinear 
Schrodinger equation describes the electromagnetic field envelope in a single polarization case. Actual optical fibers 
however, have the property that there is a difference in the propagation velocity for the two different polarization 
states of the electromagnetic field. This property is called birefringence. The effects of birefringence are almost never 
steady, but they vary randomly, in both magnitude and orientation. The random birefringence will ensure that an 
initially localized pulse will eventually disintegrate. This effect is called Polarization Mode Dispersion (PMD). PMD 
is important in situations where high-bit-rate optical signals have to be transported over long distances such as sub- 
marine intercontinental optical connections. This is the reason that the topic has already had considerable attention 
in both theoretical and experimental research in the field of optical telecommunications [0-|f3 



The dynamics of nonlinear waves in birefringent optical fibers is described in the literature by the following set of 
coupled nonlinear differential equations, which were originally introduced by Berkhoer and Zakharov [n2[: 



iU2x - i5u2t + -jU2tt = -(^2^ + l\ui\^)u2- 



(1) 



A detailed derivation of Eq.(|l|), which are from now on called the birefringence equations, can be found in the book 
by Hasegawa [14 . In Eq.([^), the (scaled) electromagnetic field envelopes of the different polarization states (modes) 
are described by ui and U2 respectively. The parameter 7 describes the strength of the cross-phase modulation. For 
single mode optical fibers we can use 7 = §• The parameter 5 describes the group velocity birefringence. In the limit 
(5 — > 0, the well known Manakov equation is retained. Ueda and Kath have argued that the parameter 5 is actually a 
random parameter with zero mean 

{6) = 0. (2) 

As a result of this, Eq.(|^) is a stochastic differential equation. In the following, we treat the parameter 5 as the 
half-dilference of the group velocity between the local principal birefringence axes. If we are able to solve Eq.(|l|) 
analytically however, we can obtain analytical expressions for the propagation properties of an initially localized 
soliton in a birefringent medium. 

The set of coupled equations ([^) have also been investigated by others. In Ref.|15| the behavior of solutions of 
Eq. (H) is investigated numerically. One of the conclusions of Ref. is that for solitonic initial conditions the two 
partial pulses (modes) always move together. A complicated pattern of soliton-soliton interactions between the modes 
is presented. Another important result is presented in Ref. . Firstly, it is shown in this publication that by applying 
the transformation: 



the set of equations Eq.(hl) transform into: 



iU2x 



ui = ui exp 
U2 = U2 exp 



1 

77"ltt 



-U2tt 



I X — lot 

2 
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\U2? + 1\U1\^)U2- 



(3) 



(4) 
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If 7 = 1, the set of equations (Q) are called the Manakov equations. The Manakov equations are well-studied and 
solutions can be obtained in several ways |l^, |9[ |2^, 22 ., 2^. It has to be remarked that the transformation 
leads to special solutions of Eq.(Q), in which the effect of the birefringence appears in the phase factor. Since under 
experimental conditions usually the optical power is measured {i.e. the transformation (|^) leads to special 

solutions of Eq.(^, which can not explain the decay of optical pulses due to birefringence (since the birefringence is in 
the phase). On the other hand in Ref.||l^ solutions are reported which show that Eq.(U) has modes which propagate 
with a different velocity. A general method to solve the equations (0) has never been given [l^ . 

There are more approaches possible which can be used to find solutions of Eq.(|l|). Although Hirota's method has 
been proved to be successful for many types of equations, we approach the problem in this paper using a method 
developed by one of us (see Ref . |Q ) to find special solutions of Eq. (0) . These solutions have the property that the 
solitonic part of the modes propagate with a different velocity. The existence of such a solution might be important 
for explaining the decay of solitonic pulses in birefringent optical fibers. For linear birefringence it is well known that 
the decay of initially localized pulses can be explained by the fact that both modes propagate with a different velocity 
along the birefringence axes [ psf . The results of this paper suggests that for solitonic pulses a similar behavior takes 
place. 

In Sect. II, we apply concepts derived in earlier publications to find special solutions of the birefringence problem. 
We derive solutions consisting of of two unperturbed solitons propagating with a different velocity along the principal 
birefringence axes. Moreover the solutions contain interactions terms between the pure solitonic solutions. It is 
concluded that two effects play a role in the decay of solitonic pulses. The effect is similar as in linear birefringence 
and is due to the unequal propagation velocity of the modes. The second effect, which is absent in linear birefringence 
is the nonlinear soliton-soliton interaction between the modes. The paper is concluded with a discussion. 



II. Special solutions of the birefringence equations 

For solving Eq.(^, we try a starting point the following free wave solutions: 

ui{x,t) = Ae*(fei*-'^i^)^ 



If we substitute Eq.(||) into Eq.(^, we obtain: 



[lUi — Ski — ^^l] yle*'-'^^*"'^^^-' = — y4|^|2g2i(fcit-(.Ji3:)g-i(feit-wi2;) _ ^j^^j^^2 ^i{kit-i.Jix) ^i{[k2--k2]t-[w2-(J2]x) 



[lU2 + 6k2 - \kl] BeK'=2«-'^2^) = _5|5|2g2j(fe2t-^2a;)g-«(fe2t-a;2a:) „ ^ Q\^J^\2 ^i(k2t~U2x) ^ 



2 i(k2t~u]2x) i([ki-ki]t-[uii-uji\x) 



(6) 



In Eq.(^, it is assumed that ki and uji are both complex (the bars indicate the complex conjugates). We require that 
the following dispersion relationships are valid: 



\ (7) 
UJ2 = -Sk2 + -k^. 

We can conclude that for nonzero coefficients A and B, the try-solution (||) does not satisfy Eq.(|^) and hence the bire- 
fringence equations have no free wave solutions. Before deriving special solutions of Eq.(^ the following simplifications 
are imposed to simplify the computations which will follow: 

/ : ki — k2 — k; k = a + ib, 
II : a — 0: k = ib. 
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Condition I is related to the fact that we strive for a situation in which both the hnear part as the nonhnear part of 
Eq.(0) are expanded in the same basis functions. In experiments, this situation can be obtained by optimizing the 
"launching conditions" . Condition II is introduced by realizing that in fiber optics only envelope solutions can be 
measured. In this paper, we search for modes ui and U2 propagating with unequal velocity. We therefore substitute 
the following solutions into Eq.(0): 



ui (x, t) = e 2 ^ ^ ian+i.me" 



[2n+l]hzi -2m6z2 



n—0 m— 

C30 C30 



(9) 



n— m— 



In Eq.(^) it is used that zi = t — Sx and Z2 = t + Sx. For the solutions (^) the left-hand side of Eq.(|l]) is equal to: 



1 



CXD OC r 



n—0 m— 

^ OQ OC 

lU2x - '^0U2t + -U2tt = 



EE 

n—0 m— ^ 



-ifc^ - AmibS + i {2n + 2m + 1)^ b'^ 
+ 4TOiM + ^ (2n + 2m + 1)^ b'^ 



; -[2n+l]bzi -2m62:2 

^2n+l,mt; o , 



(10) 



-[2n+l]6z2g-2?rifczi 



We can conclude that the linear part of Eq. is not modifying the structure of the exponential functions . This 
appears also to be the case for the nonhnear part. If we substitute Eq.(^) into the right-hand side of Eq.(|^) we obtain: 



OC OC n — 1 p m v 



/ ^ / ^ / ^ / ^ / ^ ^2g+l,u;^2p-2g+l,-u-t«^2n--2p-l,m-t)e 6 



n— 1 m— p— g— d— -u;— 
OC OC n pm — li? 



1 ^,^^'6^ R R A ^^[2n+l]f)z, -2m6z2 

+ 7^^ 2^ 2^ 2^ 2^ 2^ 2^ i'2m-2t)-l,n-pi'2u)-l,q^2p-2q+l,i;-u)e 6 



n— 1 m— 1 p— g— d— w—Q 

OC OC 1 p m V 

\ul\u2 + ^\U\\U2 = 62*'' E E E EE E ^2q+l,w>S2p — 2q-\-l^v — 'W-^2n — 2p—l.m — v^ 

n—0 m— p— g— i)— iti— 
oo OC n p m— 1 v 



a; 



EV^ 1 4 D ^-[2n+l]622^-2mfc^i 

/ ^ 2-^2-^ 2-^ 2^ ^2m-2v-l,n-p^2w-l,q-D2p~2q^l,v~w^ ' ' 6 



n—lm—lp—Oq—O v—0 w—0 



(11) 



It can be concluded by comparing Eq.(lO) and Eq.(ll) that solutions of the type (|9|) ensure that both the linear part 
and the nonlinear part of Eq.(|l]) can be expanded in the same basis functions. Given non-zero coefficients Ai.o and 
Bifi {A,B £ Ft), all the other coefficients A2n+i,m and i?2n+i,m of the solution are determined by the following 
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iteration series: 

n— 1 p m V 



— 2q-[-l.i} — W'^2n — 2p—l^m — v 



62 - -imibd + i (2n + 2m + 1) fe2 

S2m-2D-l,n-p-B2tu-l,gj42p-2g+l,D-iD 

(12) 



n p m— 1 

U U ^0 ^0 -5^' - 4"^^W + i (2n + 2m + 1)^ 62 

n—lp m V 



S2n+l,m - - XI 51 51 XI ^ 



B2q+l,wB2p-2q+l,v-wB2n~2p-l,r. 



62 + 4mi6(5 + i (2n + 2m + 1)^ 62 



p=0 9=0 u=OiD=0 2 ' ^■'--"^ < 2 

n p rn-l « ^ 4 R 

-^2m-2v-l,n-p^2w-l,q-D2p-2q+l,v-w 



U U ^0 ^0 + 4^6,5 + i (2n + 2m + 1)^ 62 ' 

where in the second terms on the right-hand side m > 0. 

The solution (^) ol Eq. is purely formal and does not give adequate insight of the behavior of solitonic solutions 
in birefringent media. In order to obtain this insight, it is convenient to separate the solutions ui(x^t) and U2{x,t) 
into two parts: 



oo oo 

^, (-y i\ — ehb^^X^ A „-[2n+l]bzi _j_ „hib'^ x ST^ ST^ A -[2n+l]bzi -2mbz2 

n—0 n—Om—1 

00 oo oo 

a, A — c^'i'b^^'Sr^ R -[2n+l]bz2 _i_ ^hb'^xST^ f> -[2n+l]bz2 ^-2m.bzi 

n—0 n=0 m—1 

In this special case, we find that given the coefficients ^2n+i.o and -B2n+i,o are given by: 

n—l P I I J 

1 _ \^ ^2q+1.0^2p-2q+1.0^2n-2p-1.0 

P R R R 

A _ -D2g+l,0i'2p-2g+l,0-t)2n-2p-1.0 



(13) 



(14) 



We firstly analyze the first part the right-hand side of Eq.(|13|). By assuming that Ai^ = A and Bi q = B all the other 
coefficients ^2^+1^0 and B2n+i.Q follow from Eq.(p^. This implies that Eq.(p^ can be reformulated in the following 
form M: 



i,{x,t) = -Aei^'"^e'^°secHb[t~5x]+^o)+e^^''''Y.T. A2n+i.n.e-^^"+'^''^'e-^^''^\ ^0 = log ( ) ; 



1 . ^A 



n—0 m—1 

00 00 



?2 



U2ix, t) ^ -Be5"''-e^«»sech(6[i + 5x] + Co) + e^*"'* X X ^2„+i,™e-[2«+i]^-^e-2™''-S 6 = - 2 log ( ^ 

n=Om=l ^ 

(15) 

The first term on the right-hand side of Eq.(^) represents unperturbed solitons. The propagation velocity of these 
solitons is determined by the birefringence coefficients 5. It follows from Eq.(^) that the two solitons ui{x,t) and 
U2{x, t) propagate with a relative velocity 2b5. This makes that even in an "ideal situation" , in which the non-solitonic 
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terms in Eq.(|l5|) are negligible an initially localized pulse is still unstable due do the unequal propagation velocity of 
the solitonic solutions along the principal birefringence axes. 

The second effect, which leads to instabilities, is the interaction between the unperturbed solitonic solutions. The 
asymptotic expansion of these instabilities is given by the sum of terms on the right-hand side of Eq.(^|). In principle, 
by using the recursion relationship jl^), we can compute all the expansion coefficients A2n+i.m and -B2n+i.m- Only in 
special cases it is possible to carry out the summation explicitly to obtain an analytical expression for the interaction. 

In Appendix A it is indicated that the summation in Eq. ([l5|) converges ii S < b. The interaction terms are 
proportional to a factor depending on b, S and 7. From this result the conclusion can be drawn that the interaction 
terms are small if (Appendix A): 



&(l + 7) 



26 -6b 



< 1. 



(16) 



If the condition ( |l^ ) is satisfied, the interaction between the modes is small compared to the solitonic part, and hence, 
the only that causes initially localized pulses to disintegrate is the unequal propagation velocity of the modes. 



III. Discussion 

We have shown that the equations (|l|) have special solutions which consist of solitonic pulses each propagating with 
a different velocity along the principal birefringence axes. The purely solitonic pulses are subject to soliton-soliton 
interactions between the modes. This indicates that there are two mechanisms in birefringent optical fibers which can 
be responsible for the decay of optical solitons in birefringent optical fibers. 

The first mechanism is exactly equal as for linear optical pulses. The solitonic optical pulses can decay because 
of a different propagation velocity of the two modes along the principal birefringence axes. It is no surprise that also 
solitonic pulses are subject to linear birefringence since the propagation velocity of the solitonic part of the solutions 
(^5|) is determined by the linear part (dispersion relations) of the differential equation (|l|). The nonlinearity in the 
differential equation (|l|) has only influence on the shape of the solutions and keeps the propagation velocity unaffected. 
If we consider for instance the propagation properties of solitons in birefringent optical fibers, we can conclude that 
for nonlinear birefringence similar criteria as for linear birefringence apply. If we remember that Eq.(|l|) is defined 
in a coordinate frame moving with velocity vq, we find that solutions ui{x,t) and U2{x,t) propagate with velocities 
vo — i:^)'^o and vq + {:^)vo. In Ref. [l3| it is shown by using a mathematical analysis based upon Ref.|^ that in this 
case the PMD in long birefringent fibers increases as the square root of the length. 

The second effect, which can lead to decay of solitonic optical pulses is the nonlinear soliton-soliton interaction 
between the modes. This effect has been reported in Refs.||l5|, In the numerical results presented in Ref.|^ it 
is concluded that only soliton-soliton interactions played a role in the decay of solitons. It is not observed that two 
initially solitonic solutions propagate with a different velocity along the birefringence axes. We think that this result 
can be explained by the fact that in Ref.|jl^ perfect solitonic pulse-shapes are taken as the initial condition. If we 
solve Eq.(^), it follows from Eq.(^ that in the case that (5 = 0, both modes can be chosen to have equal dispersion 
relations. Hence the solutions of Eq.(^ can be expanded in the following form: 

00 



00 



n=l ^^^^ 



n=l 



For nonzero coefficients Ai and B2 , all the other coefficients can be determined by a similar iteration series as presented 
in H^. This implies that the Manakov equations (El) have solitonic solutions. It might be possible that by choosing 
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purely solitonic initial conditions, the numerical simulations in Ref. [isl] lead to the special solutions (^. In realistic 
(experimental) optical systems it is difficult to make pure solitons and a certain amount of noise will always be present. 
The appearance of noise can cause that the modes propagate with a different velocity and that the actual situation is 
more complicated as presented in fl5| . 

We want to conclude by remarking that there is another way to understand the result obtained in Eq.(p^. We 
therefore firstly implement the following Galileian transformations: 



zi 

Z2 



Sx 
Sx 



As a result of this transformation we find that for e = 1, Eq.(0) can be reformulated as: 



lU2x 



-U2z 



+ W2 U2 



-£j\Ui\ U2. 



(18) 



(19) 



Eq.(|19|) describes a perturbed nonlinear Schrodinger equation for both the modes iii and U2- It should be remarked 
that both the modes ui and U2 have a different time-evolution due to the rescaled time. If we try to solve Eq.(|l^) by 
implementing a perturbation series: 



ui = /o + e/i + 6^2 + (^h + 
U2= 90 + £51 + £^52 + e^.g3 + ' 



(20) 



we obtain for the lowest order: 



«/o. + |/o.i.i + I/0P/0-O, 
igox + 2.90^222 + \9q?9o = 0. 



(21) 



19D is the nonlinear Schrodinger equation. The solutions for /o and 50 are therefore nonlinear Schrodinger solitons 
which are defined with respect to the transformed time zi = t — 5 and Z2 = t -\r 5 respectively. Furthermore, /o and go 
are contaminated with the higher order perturbations /„ and gn- This is exactly what we find in Eq.(^5|) where the 
asymptotic behavior of /„ and gn is presented. As expected, the functions /„ and gn have a time-evolution which is 
a combination of (exponential) functions which are defined with respect to the transformed time zi and Z2- 
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Appendix A 



Convergence of Eq. (|T5|) 



In this Appendix the convergence of the summation in Eq.(|l5|) is discussed. We estimate the coefficients ^2n+i,m and 
B2n+i,m with respect to A2„+i,o and i?2Ti+i,o- 



A 



71— 1 p m V 

EEEE^ 



n n n n-kb^ ~ ^mib6 + \ {2n + 2m + lY &2 

p—Q q—0 v—0 w—0 2 ' 2 ^ -' 

n p m-1 V ^ R A 

_l_ J^2m-2v-l,n-p02w-l,q-^2p-2q+l,v-w 

^ 2^2^ 2^ 2^ 1 ^2 _ 4„^i5j + i (2n + 2m + 1)^ 62 



p=0 q=0 v=0 w=0 2 
n—lp m V 



< (1 + 7) 



< (1 + 7) 



C2q+l.wC2p-2q+l,v-wC2n-2p-l,m-v 



P ™ " C2q+l.a,C2p-2q+l,v-wC2n-2p-l,ra-v ( - + |(2n + 2m + 1 )^ 6^ 

Z-^ Z-^ Z-^ Z-^ 



(Step 1) 
C'„,„=max{A„,„,B„,,„} 



=0 q=0 D=0 to=0 



-i62+i(2„ + 2m + 1)^62 _^ 16^252^2 



C2q+l,wC2p-2q+l,v-wC2n-2p-l,m-v {'imbS) 



n—1 p ni V 

EEEE 

p=o q=ov=Ow=o ( -i62 + 1 (2n + 2m + 1)'^ 62 j ^ 26m262J2 



< (1 + 7) 



(5„ „=max{i„,„,i3„,„,} 

P ™ ^ C2q+l,wC2p^2q+l.v~wC2n^2p-l,rn-v (-56^ + i (2n + 2m + 1)^ 62 

p=o q=ov=Ow=o / -^52 _|_ i (2n + 2m + 1) 62 j 



n—1 p m V 

EEEE 

p— Q— ij— w—Q 



C2q+l,wC2p-2q+l,v-wC2n-2p-l,m-v {4:rnb6) 



(-i62 + i(2n + 2m + l)2 62y 



('S'tep 3) 

c„_„=max{A„,„,i3„,„,} 



< 



6(1 + 7) 



2(5-66 
6(1 + 7) 



25-66 



ri-l p 

EE 

p=0 13=0 
^^271+ 1,0 



C2q+l,oC2p-2(j+l, 0^271-2^-1,0 



-i62 + i (2^ + 1)^62 



C2rL + l,0 — ^^^{^2n + l,0;B2n + l,Q} 



(Step 4) 

C'2n + l,0—m-^^{ ^271 + 1,0 7-^271 + 1,0} 



(A-1) 



The coefficients B2n+i,m can be estimated in a similar way as the coefficients -B2n+i,m- In Step 1, we replace either 
the coefficients An^m or i?™,™ by C„^rn = niax{A„^m, 5n,m}- This replacement is justified because: 



+ 7|u2|^wi < (1 + 'y)\vfv V = max{ui,U2} 



(A-2) 



In Step 2, the triangle inequality is applied to the complex denominator. In Step 3, is the quotient presented in in 
Eq.(A-l) reaches an upper-limit by putting m = in the denominator. Finally, in Step 4 it is used that: 



^62 + 1 (2n + 2m + 1)^ 62 > 4m65 



(A-3) 
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Eq.( A-3 ) is valid if S < b for all n,m > 0. Furthermore, the sum mati on over m is estimated to be smaller than m? 
times the summation over coefficients for which m = 0. From Eq.(A-S), it follows that: 



2(5-66 



(A-4) 



As a result of Eq.(A-l) we can estimate an upper limit of the double sum in Eq.(|T^. Because the coefficients 
|^2n+i.m| and |i?2n+i,m| a^c proportioual to |42n+i,o| and |i?2n+i,o|i we can estimate the summation over n in Eq.(|l5|) 
is proportional to the unperturbed soliton. The summation over m can be carried out independently, and remains 
finite because it forms a geometric series. Finally we obtain: 



ui{x,t) < i^e^*'''^e*«°sech(fe[t - Sx]+^o) + \c 
U2{x,t) = ^Be^'^^''e'^"sech{b[t + Sx] + |o) + 



6(14 


7) 


2^- 


66 


6(14 


-7) 


2(5- 


66 



c=niax{A,B} 
c=max{A,B} 



(A-5) 
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